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When a metal is subjected to strong magnetic field B nearly all measurable quantities exhibit 
oscillations periodic in 1/B. Such quantum oscillations represent a canonical probe of the defining 
aspect of a metal, its Fermi surface (FS). In this study we establish a new mechanism for quantum 
oscillations which requires only finite segments of a FS to exist. Oscillations periodic in 1/B occur 
if the FS segments are terminated by a pairing gap. Our results reconcile the recent breakthrough 
experiments showing quantum oscillations in a cuprate superconductor YBa2Cu306.5i, with a well- 
established result of many angle resolved photoemission (ARPES) studies which consistently indicate 
"Fermi arcs" - truncated segments of a Fermi surface - in the normal state of the cuprates. 



In conventional metals superconductivity can be un- 
derstood as a pairing instability of the Fermi surface^. 
Recent unambiguous identification of Shubnikov - de 
Haad^ and de Haas - van Alphen 4 oscillations in 
YBa2Cu3 06.5i (YBCO) in high magnetic fields ushered 
a new era in the field by furnishing a long awaited proof 
that a Fermi surface exists in a high-temperature super- 
conductor (SC). This important discovery, which accord- 
ing to the conventional paradigm implies a closed Fermi 
surface (FS), creates interesting new puzzles. The exis- 
tence of such a closed FS contradicts the well-established 
result of many angle resolved photoemission (ARPES) 
studies which consistently indicate "Fermi arcs" - trun- 
cated segments of a Fermi surface - in the normal state of 
cuprates^^l In this theoretical study we establish a new 
mechanism for quantum oscillations which requires only 
finite segments of the FS to exist. We present arguments 
that this new mechanism is relevant to the cuprates and 
show that it accounts for the quantum oscillations in a 
model that exhibits genuine Fermi arcs terminated by a 
pairing gap, consistent with the ARPES data. 

Quantum oscillation experiments in YBCO0^, 
when analyz ed using the conventional Onsager-Lifshitz 
picture^D^, indicate small Fermi pockets each with an 
area covering approximately 1.9% of the first Brillouin 
zone (BZ). Such small Fermi pockets do not arise natu- 
rally from the band structure calculations and their to- 
tal area is inconsistent with the nominal doping concen- 
tration. To explain these results proposals for various 
states with broken t ranslational symmetry have been put 
forward^ 1 * 12 I 13 I 14 E1 ; leading to complicated band struc- 
tures with multiple Fermi pockets. One would, how- 
ever, expect to see signatures of such a Fermi surface 
reconstruction in other experiments, most notably the 
ARPES, which is capable of mapping out Fermi surfaces 
of metals with high accuracy. Yet, extensive ARPES 
studies on various cuprate materials show no evidence 
for Fermi pockets. 

One can attempt to reconcile the quantum oscillation 
data with ARPES by postulating that Fermi pockets do 
exist but cannot be seen by ARPES for various reasons. 



Of these, the suppression of the photoemission intensity 
on the pocket's back side due to the coherence factors is 
often citecffl 



In this study we adopt a radically different point of 
view - we assume that the Fermi arcs observed in ARPES 
are real and ask if such genuine Fermi arcs can give rise 
to quantum oscillations. With one additional experimen- 
tally motivated assumption^, namely that the gap ter- 
minating the arcs is of the pairing origin, we find the 
answer to be affirmative. Our reasoning that underlies 
this conclusion is a natural extension of the conventional 
Onsager-Lifshitz picture supplemented by the analysis of 
what happens once the electron reaches the arc endpoint. 
These semiclassical considerations are then supplemented 
and confirmed by exact numerical diagonalizations of a 
fully quantum lattice model. 



In what follows we advance two principal ideas. First, 
we formulate a simple, experimentally motivated model 
for a normal state of underdoped cuprate superconduc- 
tors that exhibits genuine Fermi arcs terminated by a 
pairing gap. Second, we demonstrate that in an ap- 
plied external magnetic field the low-energy density of 
states (DOS) in this model oscillates as a function of 
energy as well as magnetic field. The DOS is peri- 
odic in energy with frequency that is linear in 1/B and 
proportional to the Fermi arc length. At the Fermi 
level (or any other fixed low energy) the DOS oscillates 
with 1/B. The frequency of oscillations is not related 
to any area in momentum space; instead it is propor- 
tional to the gap amplitude and the fraction of the FS 
that is gapped. The origin of the quantum oscillations 
in our model is the periodic appearance of low energy 
Andreev-type bound states^ associated with the fermi 
arcs and is distinct from all mechanisms proposed to 
explain quantum oscillati ons in cuprates in the existing 
literature^- 13 ' 14 ' 15 ' 19 ' 20 '. 
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FIG. 1: The Fermi-arc metal. a) The gap structure of the Fermi arc metal (FAM) in the shaded areas adjacent to the 
BZ faces the gap is nonzero while the excitations in the rest of the BZ are gapless. Blue lines represent the gapless segments 
of the Fermi surface known as the 'Fermi arcs'. When a perpendicular magnetic field B is applied to the system electrons in 
the gapless regions follow semiclassical orbits which coincide with contours of constant energy ek- b) The amplitude of the 
modified d-wave gap Ak, Eq. ([2]), with r — 0.01 plotted along the line connecting points (0, it) and (7r, 0) of the BZ. c) Spectral 
intensity, A(k, uj) = — 7r _1 ImG'(k, 00), of the Fermi-arc metal at 00 = 0. Here G(k,cj) = [cu — e^ + iF — A^/(cu + e^ + iF)] -1 is the 
Green's function of Hamiltonian (pi, with additional impurity broadening represented by F — 25meV. We use the tight-binding 
parametrization for ek of Ref. 21 and Ao = 50meV, v — 0.6 and r = 0.01. d) Simplified model of FAM with circular Fermi 
surface and piecewise constant gap used in semiclassical analysis. We denote the angular extent of the arcs and gapped regions 
as O a and g = 7r/2 — a respectively. 



I. FERMI-ARC METAL 

We now provide a justification for the above claims. 
Our starting point is a simple phenomenological model 
for the 'Fermi-arc metal' (FAM) that we take to describe 
the non-superconducting state of underdoped cuprates 
with the quasiparticle excitation spectrum adiabatically 
connected to the ordinary <i-wave superconductor. The 
latter is well-known to describe the superconducting 
ground state of cuprates. It has a pairing gap Ak = AoXk 
with Xk = ^( cos k x — cos k y ) and quasiparticle excitation 
spectrum 



(i) 



with ek the band dispersion referenced to the Fermi level 
ep- We assume that the excitations of the Fermi-arc 
metal have the same BCS form as Eq. ([I]) but with Ak re- 
placed by a modified d-wave gap Ak which vanishes along 
the BZ diagonals as illustrated in Fig. [IJl Although the 
detailed form of Ak is unimportant for the subsequent 
discussion we often find it useful to parameterize it as 



Xk 



o(v 2 -Xl)/T 2 



(2) 



Here v sets the length of the arc while r controls the 
width of the step between the large-gap and zero-gap 



regions as illustrated in Fig. [T}d. The ordinary d-wave 
superconductor is recovered in the limit v,t — > 0. 

Several remarks are in order before we address quan- 
tum oscillations in this model, (i) We think of the above 
postulated Fermi-arc metal spectrum as originating from 
a BCS-type pairing Hamiltonian 



+ h.c.) (3) 



ka 



where c£ creates an electron with momentum k and spin 
a. Although the microscopic Hamiltonian that would 
stabilize the mean- field state described by Eqs. ( 2|3 ) is 
presently not known, we see no fundamental reason why 
such a state could not occur for suitably chosen electron- 
electron interaction. We remark that Fermi arcs have 



been argued to appear in a phase-fluctuating Sd 22 * 23 * 
and various more exotic quantum states such as the "al- 
gebraic charge liquid" 24 , (ii) A system described by 
Hamiltonian ([3| would in fact be a superconductor, al- 
beit with low superfluid density due to ungapped por- 
tions of the Fermi surface. We thus view this BCS 
Hamiltonian as an effective mean-field theory valid on 
intermediate length scales; at long length scales (com- 
pared to the magnetic length) phase fluctuations disrupt 
long-range superconducting order and render the system 
normal, (iii) The motivation for our phenomenologi- 
cal model comes primarily from experimental consider- 
ations. Indeed many experimental studies hint at the ex- 
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istence of t he residual sup erconducting order in the pseu- 
dogap stat d 25 | 26 | 27 | 28 | 29 | 30 ( A Bogoliubov-type dispersion 
of quasiparticles in cuprates above T c , observed very re- 
cently by ARPES 17 , provides further direct evidence for 
the pairing nature of the pseudogap. The spectral in- 
tensity computed from our model and displayed in Fig. 
[]J is in detailed agreement with this data. Finally, (iv) 
we note that within the space of mean-field models the 
above construction appears to be the only way to con- 
struct genuine Fermi arcs. As mentioned above, particle- 
hole instabilities without exception produce closed Fermi 
surfaces which can only terminate at BZ boundaries. 



II. SEMICLASSICAL ANALYSIS 

In the absence of a gap, the classical equations of mo- 
tion for an electron wave packet are 

e 1 
hk = — (v k xB), r = v k = -V k e k . (4) 

c a 

Eqs. Q imply that the electrons move on constant en- 
ergy contours in momentum space. The time it takes to 
complete a cycle is T = 2tt/uj c with uo c — eB/m c c the 
cyclotron frequency. 

When constructing a wave packet to describe the semi- 
classical motion in the FAM, clearly the motion on the 
arcs is governed by Eq. (4|. However, the evolution in 
momentum space drives the wave packet into the gapped 
region where its charge is no longer a good quantum num- 
ber due to electron-hole mixing induced by the pairing 
gap. We therefore choose to construct our semiclassi- 
cal wave packet out of 'bogolons' rather than electrons. 
These are Bogoliubov-de Gennes (BdG) quasiparticles of 
the underlying superconductor and can be thought of as 
coherent mixtures of electrons and holes. The bogolon 
wavefunction \I/ k = (ii k ,v k ) is an eigenstate of the BdG 
Hamiltonian 



H k = 



e k 



= e k cr 3 + A k ai, 



(5) 



where <jj are the Pauli matrices acting in the particle-hole 
space. 

The momentum k of the bogolon wavepacket continues 
to evolve according to Eq. Q since both the velocity v k 
and the charge e have opposite signs for the particle and 
the hole components. The real-space motion, however, 
is sensitive to the particle-hole mixing since the particle 
and the hole move in opposite directions. This leads to 
a modified real-space equation of motion 



Kl 2 )v k , 



(6) 



which represents the net center of mass motion of the 
bogolon wavepacket. 

The semiclassical approximation amounts to introduc- 
ing periodic time dependence in the Hamiltonian, in lieu 
of the magnetic field: TL(t + T) = 7i(t) = 7Yk(t)- Unlike 



in the gapless case where H(t) =const, here H(t) exhibits 
time dependence due to the momentum-dependent gap 
A(t) = Ak(t), while e = e k remains a constant of motion. 
To find the solution of the time-dependent Schrodinger 
equation, 



(7) 



we employ the Floquet theorem^, which is the analog 
of the familiar Bloch theorem for time-periodic Hamilto- 
nians; it states that solutions of Eq. ^ have the form 
tf(t) = e- i§t / h f E (t) with f E (t + T) = f E (t). Here 
e -iET/h an( ^ f E (T) are the eigenvalue and the eigenstate, 
respectively, of the Floquet operator 



T = T exp 



f 

Jo 



H(t)dt 



(8) 



and T represents the time-ordering operator. If we re- 
gard the two-component structure of ^(t) as a pseu- 
dospin, then the Floquet states precess about a time de- 
pendent axis ez + A(t)x. 

The quantity E has dimensions of energy and is closely 
related to the quasiparticle energy of the original time- 
independent problem. Since the Floquet equation yields 
e -iET/n anc j t — 2ir /u c \t is clear that E is defined only 
modulo hw c . This is analogous to momentum being de- 
fined only modulo reciprocal lattice vectors in the Bloch 
theory. Henceforth we refer to E as 'quasi-energy'. 

In order to obtain analytic results we simplify our 
model further. We assume a free-electron dispersion 
e k = h 2 k 2 /2m — ep and that the FAM gap is piecewise 
constant and dependent only on the momentum direc- 
tion. We take A k equal to zero on the arcs and ±A 
elsewhere as illustrated in Fig. [l]i. The details of the cal- 
culation of the Floquet quasi-energy are provided in the 
Methods. Here we present the results and discuss their 
implications. 

The quasi-energy as a function of both the band energy 
e and the magnetic field is shown in Fig.[2^i. It contains 
most of the physics of this model. Fig. [2p shows a cut 
along the e direction for constant B. The quasiparti- 
cle energy dispersion is obtained by a simple 'unwinding' 
procedure, described in Methods. Energy bands sepa- 
rated by small gaps result, in close analogy to the Bloch 
energy bands. The density of states in Fig. ^ displays 
clear periodic structure with frequency that can be esti- 
mated from Eq. (16) as 



F e = 



(9) 



This is in agreement with the exact numerical results 
which are discussed in the next section. 

We now turn our attention to the low-energy behav- 
ior of the FAM as a function of field B. Near the Fermi 
energy, e — > 0, the quasi-energy coincides with the quasi- 
particle energy, E = E, and no unwinding is necessary. 
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FIG. 2: Semiclassical analysis, a) Density plot of the quasi-energy E as a function of the inverse magnetic field (times 
Aora c c/e) and the band energy (in units of the gap amplitude). E is analytically obtained for the simplified Fermi-arc metal 
with a piecewise constant gap structure; it is scaled by hco c and ranges between (dark) and 1 (bright), b) Quasiparticle 
excitation energy E, in the Fermi- arc metal as a function of the band energy e at constant B. The quasi-particle energy is 
related to the quasi-energy E as E = E + nhu> c , where n, the unwinding index, is chosen to continuously match E with e 
for e ^> Ao, and maintain E as a piecewise continuous and monotonically increasing function of e. c) Density of states at a 
constant magnetic field B w 40T as a function of energy, d) The quasi-energy E, for several low lying states with e < hco c 
vs. 1/B. e) The density of states at the Fermi energy vs. the magnetic field. In both plots the given field range corresponds 
to 20 — 60T. The number of DOS oscillations found in this range is about 17, in good agreement with the experiments. In all 
panels a — tt/4 corresponding to v ~ 0.6 in the numerical study. 



When e — > the quasi-energy vanishes, but the den- 
sity of states, which is proportional to dE/de, depends 
strongly on the magnetic field. In Fig. [2]i we present the 
quasi-energy for a few different values of e close to the 
Fermi energy as a function of 1/B. In certain magnetic 
fields the density of lines is high and this translates to the 
sharp peaks in the DOS shown in Fig. [2^ . This result 
is directly related to the experimentally observed oscilla- 
tions. From Eq. (16) we may deduce that a peak in the 
Fermi energy DOS occurs whenever cos(Ao# g /hu; c ) = 0, 
leading to oscillation frequency 



F = An 



it he 



(10) 



Using Aq = 80meV, g = 7r/4 (Ref. [32]) and a cyclotron 
mass m c = 3m e (Ref. 3) we estimate F = 518T, very 
close to the dominant experimental frequency 530-540T 
in YBCO. 

Some intuitive understanding of the origin of the DOS 
oscillations can be gained by examining a typical low- 
energy Floquet state and its associated real-space 
trajectory. This is illustrated in Fig. [3] The states con- 
tributing to high DOS are reminiscent of the Andreev 
bound states found on extended impurities and on sam- 
ple edges in d-wave superconductors^. The periodicity 
in inverse magnetic field in this model is a consequence 
of the periodic appearance of these Andreev-type states 
on the Fermi arcs at low energies. 

The quantum oscillation mechanism described above 
is dominant for quasiparticle energies much smaller than 
the gap amplitude Aq, and does not involve the conven- 
tional Onsager-Lifshitz action considerations. Obviously, 



the action must be considered to obtain the normal-metal 
behavior in the limit of small gap or large arcs. The in- 
terplay of the two mechanisms then becomes complicated 
and we briefly discuss it in Methods. 



III. LATTICE MODEL 

In order to exemplify the validity of our semiclassi- 
cal analysis we now consider a fully quantum-mechanical 
lattice model of the Fermi-arc metal and confirm the ex- 
istence of quantum oscillations by exact numerical calcu- 
lations. To this end, we study the real-space version of 
the Hamiltonian (J3|, 



w = 5>. 



j A c 



i j a 



(ii) 

where c\ a creates an electron with spin a at site of the 
square lattice. The effect of magnetic field is described by 
the usual Peierls factors Oij = (2tt/^o) J r J A • dl and the 

SC order parameter A(r^, r^) is correspondingly taken to 
contain a periodic lattice of Abrikosov vortices. This is 
achieved by adopting 



A(n,r 2 ) = e 



,B(R)^ e -k Ak5 

k 



(12) 



where R = (ri + r 2 )/2 is the center of mass and r = 
ri — T2 the relative coordinate of the Cooper pair. The 
vortex lattice is reflected in the phase O(R) winding by 
2tt around each vortex while the Fermi-arc structure in 
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FIG. 3: Semiclassical trajectories. Left column: Sample real-space trajectories. Right column: the quasi-particle group 
velocity relative to vf, \u\ 2 — |^| 2 , as a function of time along the trajectory. When e <C Ao the bogolon pseudospin lies in 
the y — z plane. Its projection on the z axis is \u\ 2 — \v\ 2 . The states contributing to high DOS [panels (a) and (b) with 
huj c — Ao#p/(4 + |)7r] have their pseudospin pointing in the y-direction on the arc and are precessing in the y — z plane in 
the gapped region. This means that on the arc the state is a perfect mixture of a particle and a hole, resembling a familiar 
Andreev bound state. The associated real-space trajectory encloses a particularly small area, since on the arc, the group 
velocity vanishes, \u\ 2 — \v\ 2 = 0. On the other hand, the states near the DOS minimum [panels (c) and (d), hcu c = AoOg/An], 
point in the z pseudo-spin direction and are thus either particles or holes. Their motion is dominated by the arc region, and 
the real-space trajectory encloses a large area. The DOS oscillations with B may be understood as follows. In order to exhibit 
no real-space motion on the arc, Andreev states must point in the y pseudospin direction. In the gapped region, the pseudospin 
precesses about the x-axis (assuming e« Ao). If the magnetic field is such that the time to traverse the gapped region, 9 /uj c , 
allows exactly an integer plus one half cycles, the pseudospin will end up in the — y direction and will continue to have zero 
motion on the next arc, thus giving a consistent Floquet state. 



/c-space follows from taking Ak as described by Eq. ([5]). 
In order to keep subsequent calculations simple we limit 
ourselves to the nearest-neighbor electron hopping in the 
kinetic term: tu = ep, Uj = —t for nearest neigh- 

bors and tij = otherwise. 

We solve the problem posed by Eqs. ( |11|12 ) by employ- 
ing the Franz- Tesanovic (FT) transformation 33 . This 
unitary transformation removes the non-trivial phase 
0(R) from the pairing term and renders the transformed 
Hamiltonian translationally invariant with a unit cell 
containing two superconducting vortices. Eigenstates of 
this new Hamiltonian then can be conveniently found by 
appealing to the Bloch theorem. Ref. [34] gives a detailed 
description of the implementation of the FT transfor- 
mation to the lattice model of s-, p- and d-wave super- 
conductors. The treatment of our modified d-wave SC 
follows as a straightforward generalization of this proce- 
dure. 

The results of our numerical calculations are summa- 
rized in Fig. [4] The density of states D(e) of the mod- 
ified (i-wave SC in the applied magnetic field shows the 
expected Landau level structure^ at energies exceeding 
about 3A . The surprising new result is the appearance 
of a clear periodic structure in D(e) even at low energies 
inside the SC gap. This is unlike the ordinary d-wave SC 
where the Landau level type oscillations are known to be 



absent at low energies^ 33 * 34 * 35 * 36 ! (we have confirmed this 
result by setting v = in our model). We attribute the 
low-energy oscillations to the gapless regions of the BZ 
implied by the modified d-wave order parameter Eq. ([2| 
with v 7^ 0. 

The power spectrum of the low-energy DOS displayed 
in Fig. |4j3 confirms the periodic structure with a period 
that scales with 1/5, in accordance with the semiclassical 
result Eq. Physical observables, such as the specific 
heat and resistivity, depend on DOS at the Fermi level, 
D(cf)- Fig. shows that D{ep) exhibits similar oscil- 
latory behavior. 

The above exact numerical results show unambiguous 
evidence for quantum oscillations in FAM, a system that 
by construction exhibits genuine Fermi arcs terminated 
by a pairing gap. Although the exact diagonalization 
technique does not allow us to study these oscillations as 
a function of smoothly varying field B (and thus com- 
pare directly to experiment), we have verified that oscil- 
lations in the energy variable are in qualitative as well as 
semi-quantitative agreement with the semiclassical pic- 
ture presented in Sec. II. Specifically, we have analyzed 
the dependence of oscillation frequency s on the gap am- 
plitude Ao and the arc length v and found these in agree- 
ment with the semiclassical predictions. This comparison 
is discussed more fully in the Methods. 
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FIG. 4: Exact diagonalization of the lattice model, a) DOS as a function of energy in the Fermi-arc metal in zero (red) 
and non-zero (black) magnetic field corresponding to two vortices in a 20 x 20 magnetic unit cell. In YBCO with lattice constant 
do — 4A this corresponds to the physical field of about 64T. The parameters used are as follows, Ao/t = 1, = —1.3, v — 0.6 
and t = 0.1. b) The low-energy DOS for the same parameters, in detail, c) The power spectrum of the low-energy DOS showing 
dominant frequency of oscillations 210£ -1 and its second harmonic, d) DOS at the Fermi level as a function of cf- 



IV. OUTLOOK 



The original observation of quantum oscillations in 
YBCO 2 has been interpreted as quantitatively consis- 
tent with the ARPES measurements^ by assuming that 
the Fermi arc observed in ARPES was a part of a Fermi 
pocket resulting from the Fermi surface reconstruction 
due to a symmetry breaking instability with wavevector 
(7r, 7r). Such a Fermi pocket would have an area of about 
2% of the BZ, consistent with the observed oscillation 
frequency. The ARPES data used for this comparison 
however pertain not to YBCO but to a different high-T c 
compound Na2- :c Ca :c Cu202Cl2 (NaCCOC). If one uses 
instead the ARPES data on YBCO that became avail- 
able more recently 32 , then the agreement disappears: the 
corresponding Fermi pocket would comprise about 5% of 
the BZ, leading to the frequency more than twice that 
observed in experiment. The two experiments are eas- 
ily reconciled by appealing to the Fermi arc picture ad- 
vocated above, where the frequency of oscillations does 
not relate to any Fermi surface area but originates from 
the periodic appearance of Andreev-type states associ- 
ated with Fermi arcs. The oscillation frequency ( 10 ) de- 
pends on the gap amplitude and the size of the gapped 
region in the Fermi-arc metal. 

The above comparison illustrates some key differences 
between the conventional Onsager-Lifshitz picture of 
quantum oscillations in YBCO— * * 13 * and the mechanism 
invoking genuine Fermi arcs terminated by a pairing gap 
proposed in this study. The latter relies only on the FS 
structure that is directly seen in ARPES while the for- 
mer must make assumptions about unseen portions of 
the FS in the pa rts of t he BZ where the ARPES indi- 
cates a large gap) 11 * 12 * ^ . A direct observation of Fermi 
arcs at low temperatures in high magnetic fields would 



discriminate between the two pictures. Under such condi- 
tions ARPES experiments are not feasible but it should 
be possible to image the FS by means of the scanning 
tunneling probe using th e techni que of Fourier-transform 
interference spectroscop y 37 * 38 * 3 -^. Even if confirmed, the 
microscopic origin of the Fermi arc phenomenon remains 
an open question, the answer to which may pave the road 
towards the full solution of the cuprate mystery. 
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VI. METHODS 

A. Experimental considerations 

A crucial assumption we made in de ducing the oscil- 
latory behavior of D(e) in Eqs. ( |9|1Q[ ) is that the gap 
structure itself is field-independent. This is equivalent 
to assuming that magnetic field enters the underlying 
mean-field Hamiltonian via the usual minimum substi- 
tution but does not alter any of its parameters. This 
known to be true in ordinary metals but it is less obvi- 
ous that this assumption applies to our Fermi-arc metal 
which relies on the existence of the residual SC gap. The 
amplitude and the /c-space structure of the latter could be 
susceptible to magnetic field. Specifically, if the Fermi arc 
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length depended on B then the frequency of the quantum 
oscillations would itself become a function of £?, in con- 
tradiction to experimental finding of constant frequency. 
Since there exist no independent measures of the Fermi 
arc length as a function of B, and since the microscopic 
theory underlying the Fermi arc phenomenon is also un- 
known, we must adopt the arc length independence on 
B in the experimentally relevant interval 20-60T as an 
additional assumption of our phenomenological model. 



B. Floquet eigenstates 

With a piecewise constant gap A(t) and a constant 
band energy e the time ordered exponent of the Floquet 
operator ^ may be written as a product of 8 operators, 
T = U-UqU^UqU-UqU^Uq. These correspond to the 
time evolution operators on the eight segments of the 
Fermi surface (4 arcs and 4 gapped segments), 



U 



-iea 3 T /h 



u± 



-i(e<7 3 ±A <7i)Ti/fi 



(13) 



where T = a /uo c and T\ = g /uo c are the times to tra- 
verse the arc and the gapped regions, respectively. The 
properties of the Pauli matrices allow us to write 



U 

U± 



COS 



( Ek@g 
cos — — - 

V nw c 



i sin 



fajJr 



i sin — — l 



(73 (14) 
e(T 3 ± A <7i 



E k 



Aq. Since the motion on the first 4 



with 

segments is the same as on the last 4 segments we may 
define T\j 2 = T and find the eigenstates of half the cycle. 
Multiplying the four operators yields 



A 

iB* 



iB 

A* 



(15) 



with 



A 



A = | — smai | + e 



COS OL\ — i —— sm OL\ 

E^ 



B 



e ia ° sm ft] ( sm a\ cos ao + cos a\ sm ao ) , 



where ao = eO a /hco c and ql\ = E\^O g /hx} c . Using the fact 
that T\/2 is unitary (and therefore |A| 2 + |£>| 2 = 1) we 
find the real part of its eigenvalue as cosine of the phase 
(j) = (27rE/hco c ) and deduce the quasi-energy 



E — arccos [Re (A)] 



(16) 



This is presented in Fig. [2k. 



C. Energy unwinding 

The full quasiparticle energy, E, is deduced from the 
quasienergy as a function of the band energy, E(e), 



through the following unwinding procedure. First, since 
the inverse cosine function in the quasienergy gives an- 
gles between and ir we interpret decreasing segments 
of the quasienergy as resulting from angles between tt 
and 27r. On these segments we replace: arccos [Re( A)] — > 
2tt — arccos [Re (A)]. The result are disconnected mono- 
tonically increasing segments of energy between and 
hw Cl the energy bands. In analogy with Bloch bands, 
each band begins and ends with a flat dispersion. We are 
free to shift the energy by an integer multiple of hw c and 
do so uniformly in each band, i.e., without 'breaking' the 
bands. Different bands are shifted by different amounts 
in order to create a monotonically increasing function. 
We define E = E + nhw c where n is the band unwinding 
index. In the high energy region, e ^> Aq, the SC gap 
is negligible and we expect the energy to converge to the 
band energy. Indeed, if we choose the band unwinding 
index to simply count the bands starting at zero for the 
first band, we obtain a monotonic function with small 
gaps separating the bands at low e. At high e we recover 
the desired linear dependence with slope 1. However, a 
small offset between E(e) and e appears. We interpret 
this offset as requiring larger gaps at low energy. Usu- 
ally steps of 2hw c or 3Huj c are sufficient. We note that 
with the exception of vanishing arc length, the slope at 
e = is finite, meaning that the Fermi energy is a band 
midpoint where a gap does not open. In addition, we 
expect the energy to be antisymmetric with respect to 
e and therefore E(e = 0) = 0. As a result, the lowest 
energy band is not shifted. 



D. Action considerations 

The quantum oscillations described in Sec. II do not 
originate from and are not affected by the action quanti- 
zation. Nevertheless, the action plays an important role 
in the limit of small gap or large arcs, when Andreev 
bound states are rare. We use the circular momentum- 
space trajectory and the Floquet state pseudospin to 
determine the real-space paths and to calculate the as- 
sociated action. At energies close to ep the action is 
proportional to the energy. To a good approximation 
S/H ex lgAk\ cos( Ao9 g /huj c )\ wrtn the magnetic length 
Ib and Ak the area in /c-space bounded by the trajec- 
tory. Thus, the linear dependence of the action on 1/B 
is modulated by | cos(A 6^/fkj c )|, the projection of the 
pseudospin on the z-direction while on the arc. There- 
fore, for magnetic fields with peaked DOS the action is 
zero. This means that the Andreev states are allowed by 
the action quantization. The action slope close to these 
points is very large (about 3.5 larger than the slope as- 
sociated with the full Fermi surface area) so that many 
other states at near magnetic fields are allowed by the ac- 
tion quantization. The fact that more states are allowed 
does not change the DOS periodicity in 1/5, however 
it may broaden the observed oscillation frequency range. 
When the gapped region shrinks (arc length increases) or 
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FIG. 5: Oscillation frequency in the lattice model. The frequency s is extracted from the power spectrum analysis 
of oscillations in the low-energy quasiparticle DOS as a function of the arc length v with a fixed chemical potential. Panel 
(a) shows scaling with magnetic field, obtained by considering 16 x 16 and 20 x 20 magnetic unit cell. Panel (b) shows the 
dependence on the gap amplitude at constant field B. 



the gap amplitude decreases the frequency of DOS mod- 
ulations due to the Andreev states will decrease. The ac- 
tion quantization then becomes more important. A sec- 
ond frequency which reflects the full Fermi surface area 
will appear and eventually dominate when the gap closes 
(the cosine is simply 1 in this limit). In this way our 
model recovers the usual Onsager-Lifshitz quantization 
in the limit of vanishing gap. 



E. Frequency analysis 

We have analyzed the behavior of the oscillation fre- 
quency s in the lattice model of Sec. Ill as a function of 
magnetic field B, arc length v and maximum gap size Ao 
(at a fixed chemical potential). Our results are summa- 



rized in Fig. [5] The oscillation frequency scales linearly 
with 1/B as expected on the basis of semiclassical Eq. 
([9]). For small and intermediate arc length v it increases 
linearly with also in accord with Eq. ([9|. Panel (b) 
shows that the frequency also exhibits dependence of the 
gap amplitude, not expected on the basis of Eq. ([9|. This 
dependence, however, is seen to saturate in the limit of 
large Ao. This is consistent with the fact that Eq. ([9| 
is valid only in the regime e <C Aq and this condition is 
not fully met in the lattice model for small Ao/t since 
we must consider a finite energy window to extract the 
oscillation frequency. A more careful analysis of Eq. 16 
indeed yields a weak gap dependence of the oscillation 
frequency F e with a trend resembling that displayed in 
Fig.|b. 
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